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Chandra Observations of Low Mass X-ray Binaries and Diffuse 
Gas in the Early- Type Galaxies NGC 4365 and NGC 4382 (M85) 

Gregory R. Sivakoff 1 , Craig L. Sarazin 1 , and Jimmy A. Irwin 2 

ABSTRACT 

We used the Chandra X-ray Observatory ACIS S3 to image the X-ray faint 
elliptical galaxy NGC 4365 and lenticular galaxy NGC 4382. The observations 
resolve much of the X-ray emission into 99 and 58 sources, respectively, most of 
which are low-mass X-ray binaries (LMXBs) associated with each of the galaxies. 
Within one effective radius of NGC 4365, about 45% of the counts are resolved 
into sources, 30% are attributed to unresolved LMXBs, and 25% are attributed 
to diffuse gas. Within two effective radii of NGC 4382, about 22% of the counts 
are resolved into sources, 33% are attributed to unresolved LMXBs, and 45% 
are attributed to diffuse gas. We identify 18 out of the 37 X-ray sources in a 
central field in NGC 4365 with globular clusters. The luminosity functions of the 
resolved sources for both galaxies are best fit with cutoff power-laws whose cutoff 
luminosity is ~ 0.9 — 3.1 x 10 39 ergs s -1 . These luminosities are much larger than 
those previously measured for similar galaxies; we do not find evidence for a break 
in the luminosity function at the Eddington luminosity of a 1.4 M & neutron star. 
The spatial distributions of the resolved sources for both galaxies are broader than 
the distribution of optical stars. In both galaxies, a hard power-law model fits 
the summed spectrum of all of the sources. The unresolved emission is best fit by 
the sum of a soft mekal model representing emission from diffuse gas, and a hard 
power-law, presumed to be from unresolved LMXBs. There is some evidence that 
the temperature of the diffuse gas increases with increasing radius. A standard 
beta model fits the radial distribution of the diffuse gas in both galaxies. In the 
elliptical NGC 4365, the best-fit core radius is very small, while the SO galaxy 
NGC 4382 has a larger core radius. This may indicate that the gas in NGC 4382 
is rotating significantly. 
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1. Introduction 

Early-type galaxies are luminous X-ray sources (Forman, Jones, & Tucker 1985). For 
galaxies of a given optical luminosity, the X-ray-to-optical luminosity ratio (L X /L B ) ranges 
over two orders of magnitude for galaxies (Canizares, Fabbiano, & Trinchieri 1987; White & 
Davis 1997). We will refer to galaxies with relatively high Lx/Lb ratios as "X-ray bright" 
and to galaxies with relatively low Lx/Lb ratios as "X-ray faint." Hot (kT ~ 1 keV) 
interstellar gas dominates the X-ray emission in X-ray bright galaxies (e.g., Forman et al. 
1985; Trinchieri, Fabbiano, & Canizares 1986); whereas X-ray faint galaxies exhibit two 
distinct spectral components: a hard (~5-10 keV) component (Matsumoto et al. 1997) and 
a very soft (~ 0.2 keV) component (Fabbiano, Kim, & Trinchieri 1994; Pellegrini 1994; Kim 
et al. 1996). Since the hard component is actually found in both X-ray bright and X-ray 
faint early-type galaxies, with strengths roughly proportional to the optical luminosity of 
the galaxy, Kim et al. (1992) suggested that the hard component is due to low-mass X-ray 
binaries (LMXBs) like those observed in the Milky Way's bulge. 

Chandra observations of the X-ray faint elliptical, NGC 4697 (Sarazin, Irwin, & Bregman 
2000, 2001), resolved the majority of emission into X-ray point sources, whose properties 
are consistent with LMXBs. Chandra observations show that a significant fraction of the 
LMXBs are coincident with globular clusters (Sarazin et al. 2000; Angelini et al. 2001). 
Taken together, the X-ray spectra of LMXBs in early-type galaxies are hard, and can be fit 
by thermal bremsstrahlung with a temperature of ~ 7 keV. However, there is a considerable 
variety in observed X-ray spectra and colors of individual sources. Some of the LMXBs 
have very soft X-ray spectra, similar to those of Galactic supersoft sources (Sarazin et al. 

2000) . The luminosity function of LMXBs in early- type galaxies in the luminosity range 
5 x 10 37 — 10 39 erg s" 1 appears to generally be fit by a broken power-law (Sarazin et al. 2000, 

2001) , with a break luminosity which is comparable to the Eddington luminosity for spherical 
accretion onto a 1.4 M & neutron star. This suggests that the sources more luminosity than 
~ 2 x 10 38 erg s -1 are accreting black holes. 

In addition to the LMXBs, the Chandra observations of NGC 4697 also detected diffuse 
emission from interstellar gas. The gas is quite cool, with temperatures ~ 0.3 keV (Sarazin 
et al. 2000, 2001). 

With Chandra's ability to study both the LMXB population and diffuse X-ray emission 
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of early-type galaxies, it is useful to extend such studies to other early-type galaxies. Since 
diffuse emission could dominate the LMXBs in X-ray bright galaxies, the ideal targets for 
studying the LMXB population of early-type galaxies are X-ray faint galaxies. In this paper, 
we present the results of two such observations. Both NGC 4365 and NGC 4382 are early- 
type galaxies located in the Virgo cluster. NGC 4365 is an E3 galaxy. NGC 4382 (M85) 
is a SO galaxy with an interacting companion NGC 4394, which together with several other 
nearby galaxies (VCC 797, IC 3292) form a galaxy group RSCG 54 (Barton et al. 1996). 
Since NGC 4382 has somewhat bluer colors than are typical for a SO, and its disk shows a 
spiral pattern, it is sometimes classified as an Sa. It was the host of the Type la supernova 
SN1960r (e.g., Porter 1993). The distances to these galaxies are 20.4 Mpc (NGC 4365) and 
18.5 Mpc (NGC 4382), based on the method of surface brightness fluctuations (Tonry et al. 
2001). For comparison, NGC 4697 has a distance of 11.7 Mpc using the same technique. 
All three galaxies have similar X-ray luminosities and L X /L B ratios (O'Sullivan, Forbes, & 
Ponman 2001). 

In § 2, we discuss the observations and data reduction of these galaxies. The X-ray 
images are presented in § 3. The properties of resolved sources are given and discussed in 
§ 4. We discuss the X-ray spectral properties in § 5. The spatial distribution of the diffuse 
X-ray emitting gas is determined in § 6. We summarize our conclusions in § 7. 

2. Observation and Data Reduction 

NGC 4365 was observed on 2001 June 2-3 with a live exposure of 40,429 s, while 
NGC 4382 was observed on 2001 May 29-30 for 39,749 s. Both galaxies were observed with 
the ACIS-235678 chips operated at a temperature of —120 C and with a frame time of 3.2 s. 
We determined the pointings so that the entire galaxy was located on the S3 chip and so 
that the center of each galaxy was not on a node boundary of the chip. Although a number 
of serendipitous sources were seen on the other chips, the analysis of both galaxies in this 
paper will be based on data from the S3 chip alone. The data were telemetered in Faint 
mode, and only events with ASCA grades of 0,2,3,4, and 6 were included. Photon energies 
were determined using the gain file acisD2000-08-12gainN0003.fits. We excluded bad pixels, 
bad columns, and columns adjacent to bad columns or chip node boundaries. 

Chandra is known to encounter periods of high background ( "background flares" ), which 
especially affect the backside-illuminated SI and S3 chips 3 . We determined the background 
count rate for each galaxy, using the SI chip to avoid the enhanced flux due to the galaxy 
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on the S3 chip. For both NGC 4365 and NGC 4382, no background flares were seen. The 
backgrounds for extended regions were taken from the deep blank sky backgrounds compiled 
by Maxim Markevitch 3 , and were adjusted to the aspect histories of our observations using 
his MAKE_ACISBG program 3 The normalization of the background was increased very slightly 
to match the 5-10 keV count rate of the blank-sky background with a galaxy emission-free 
region of each galaxy on the S3 chip. We included the "background" due to the readout 
artifact in ACIS using his make_readout_bg program 3 . 

On the S3 image of NGC 4365, there were two X-ray sources with positions that agreed 
with positions from the U.S. Naval Observatory A2.0 (USNOA2 Monet et al. 1998) or the 
Two Micron All Sky Survey (2MASS Cutri et al. 2001) catalogs. Src. 98 in Table 1, an 
X-ray source near the edge of the field where the Chandra point-spread-function (PSF) is 
quite broad and the X-ray position is uncertain, corresponded to USNOA2 0900-07143676. 
For this source, the optical and X-ray positions agree to within ~0'/4. The second possible 
optical identification is between the optical nucleus of NGC 4365 itself and Src. 1 in Table 1. 
Unfortunately, the central region of this galaxy is crowded with sources, and it is uncertain 
whether the source closest to the nucleus is, in fact, an active galactic nucleus (AGN), a 
nearly coincident X-ray binary, or a combination of sources. Confusion limits the accuracy 
with which the positions of sources near the center of NGC 4365 can be determined. In the 
absence of more direct evidence, we will assume that the positions in NGC 4365 are accurate 
to 1". In NGC 4382, there were two X-ray sources with optical identifications and positions 
from the Tycho-2 catalog (H0g et al. 2000), the 2MASS catalog, and/or the USNOA2 optical 
catalog. The first source, Src. 38 in Table 2, corresponds to J1225+182, a z = 1.19 quasar 
(Knezek & Bregman 1998) and is the brightest X-ray source in the S3 field. The second 
source, which was just below the detection threshold for inclusion in Table 2, corresponds to 
BD+18 2609, a G5 star with a proper motion of ~90 mas yr _1 (H0g et al. 1998). The X-ray 
positions agreed with the optical positions to better than 0'/4 in both cases with no evidence 
for a systematic offset. We believe that the X-ray positions for NGC 4382 are generally 
accurate to ~0'/5. 

We performed initial data reduction and some of the data analysis using the Chandra 
analysis package CIAO 4 . We extracted spectra using the calcrmf 5 software kindly provided 
by Alexey Vikhlinin and Jonathon McDowell. The spectra were fit using xspec 6 . Since 

4 See http://asc.harvard.edu/ciao/. 

5 See http:/ /asc. harvard.edu/cont-soft/software/. 

6 See http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/. 
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launch, the ACIS quantum efficiency (QE) has undergone continuous degradation 7 . We 
used the XSPEC ACISABS model to correct X-ray spectra for this effect. 

3. X-ray Image 

The X-ray images of NGC 4365 and NGC 4382 show both resolved point sources and 
diffuse emission. To best image this combination, we adaptively smoothed the Chandra 
S3 image using a minimum signal-to-noise ratio (SNR) per smoothing beam of 3. After 
correcting for exposure and background, we removed regions of less than 20 ks to avoid 
artifacts at the chip's edge. Figure 1 displays gray scale images with a logarithmic surface 
brightness scale ranging between 1.56 x 10~ 6 cnt pix" 1 s _1 and 3.3 x 10~ 4 cnt pix" 1 s _1 . In 
both galaxies, the majority of the X-ray luminosity emission has been resolved into point 
sources. 

In Figure 2, we show the Digital Sky Survey (DSS) optical images centered on NGC 4365 
and NGC 4382 using a logarithmic gray scale. The overlaid circles indicate the positions of 
the X-ray point sources listed in Tables 1 and 2. Only a small fraction of the weaker point 
sources are not evident in Figure 1. In NGC 4365, the diffuse emission qualitatively has the 
same ellipticity and position angle (PA) as the optical emission. On the other hand, the 
X-ray emission in NGC 4382 appears rounder than the optical emission. 

4. Resolved Sources 
4.1. Detections 

The discrete X-ray source populations on the ACIS S3 images of the two galaxies were 
determined using using the wavelet detection algorithm (ciao wavdetect 5 program), and 
were confirmed with a local cell detection method and by visual inspection. The wavelet 
source detection threshold was set at 10~ 6 , which implies that <1 false source (due to a 
statistical fluctuation in the background) would be detected in the entire S3 image. We 
further restricted the sources by requiring sufficient source counts to determine the source 
flux at the > 3a level. Over most of the image, the minimum detectable flux was about 
2.7 x 1CT 4 cnt s' 1 in the 0.3-10 keV band for both galaxies. The detection limit is slightly 
higher at large distances where the Point-Spread-Function (PSF) is larger, and near the 
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center of the galaxy where confusion with other discrete sources or diffuse gaseous emission 
may affect the sensitivity. For all source analyses, we used a local background with an 
area three times that of each source's region from wavdetect. There were a few cases of 
nearby sources where the source regions or background regions were altered slightly to avoid 
overlaps, without affecting the net count rates. 

Tables f and 2 list the discrete sources detected by this technique in NGC 4365 and 
NGC 4382, respectively. In each case, the sources are ordered with increasing projected 
distance d from the center of the galaxy. Columns 1-7 provide the source number, the IAU 
name, the source position (J2000), the projected distance d from the center of each galaxy, 
the count rate with its la error, and the SNR for the count rate. The fluxes were corrected 
for exposure and the instrument PSF. Since we did not detect a distinct source which was 
unambiguously associated with an AGN at the center of either NGC 4365 or NGC 4382 
(see § 4.2), we adopted the central optical/IR positions of R.A. = 12 h 24 m 28 f !26 and Dec. = 
+7°I9'3 , . / 8 for NGC 4365, and R.A. = 12 h 25 m 24!04 and Dec. = +I8°II'25'. / 9 for NGC 4382 
from the 2MASS survey (Cutri et al. 2001). As noted in § 2, the overall absolute errors are 
probably <1" for NGC 4365 and ~0"5 for NGC 4382 near the centers of the fields, with 
larger errors further out for both observations. 

Our detection limit for sources should result in <1 false source (due to a statistical 
fluctuation) in the entire S3 field of view for each galaxy. However, some of the detected 
sources may be unrelated foreground or (more likely) background objects. Based on the 
source counts in Brandt et al. (2000) and Mushotzky et al. (2000), we would expect ~12 
such unrelated sources in each observation. These should be spread out fairly uniformly 
over each S3 image (Figure 1), except for the reduced sensitivity at the outer edges of the 
field due to reduced exposure and increased PSF. Thus, the unrelated sources should mainly 
be found at larger distances from the optical centers of the galaxies (the latter entries of 
Tables 1 & 2), while the sources associated with the galaxies should be concentrated to their 
centers. 



4.2. Identifications 

In addition to the identifications mentioned in § 2, we visually examined the first and 
second generation DSS images for position matches. Both Src. 98 (USNOA2 0900-07143676) 
of NGC 4365, and Src. 38 (J1225+182) of NGC 4382 were visible. In addition, Srcs. 88 and 
97 of NGC 4365 and Srcs. 48, 53, and 57 of NGC 4382 had possible optical counterparts. 
None of these five sources corresponded to objects listed in NED or SIMBAD. 
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In both galaxies, the area enclosed by Src. 1 is coincident with our adopted optical/IR 
positions of the galaxy center. In NGC 4365, the centroid of Src. 1 is off by an amount which 
is consistent with our positional uncertainty, while in NGC 4382, the centroid is off by > 3 
times our positional uncertainty. It is unclear whether either of these sources is an AGN 
in the galaxy center or a LMXB projected near the center. We also examined hard-band 
(2-10 keV) images of both galaxies in which AGNs might be more obvious, particularly 
as the diffuse emission from the galaxies is soft. We did not see any new sources near 
the center, nor were either of the centermost sources particularly hard. Thus, we have no 
unambiguous evidence for AGNs in these galaxies. Conservatively, we adopt the luminosity 
of the centermost source in each galaxy (1.5 x 10 38 ergs s -1 in NGC 4365 and 9.0 x 10 37 ergs 
s -1 in NGC 4382) as the upper limit to the central AGN luminosity. 

We examined NGC 4382 for a source associated with SNI960r. Src. 27 is the closest 
detected source; however, it is ~20"from the Porter (1993) position. We placed a 2.0"radius 
circular aperture, similar to that of Src. 27 plus 0.5"to account for our positional uncertainty, 
and a 4.0"radius circular aperture to determine the local background. No net counts were 
detected and we can place a 3cr upper limit of 6.3 x 10 37 ergs s -1 . 

For NGC 4365, A. Kundu kindly provided a list of 325 globular cluster positions from 
data presented in Kundu & Whitmore (2001). The WFPC2 field of view overlaps 37 X-ray 
point sources from Table 1. After solving for a small offset, 18 of the 37 sources have a 
globular cluster within l'/02. The expected number of associations from random positions is 
~ 0.6 based on the density of the globular clusters. This number is consistent with Kundu 
et al. (2003), who adopted a less strict source detection criteria. Recently Larsen et al. (2003) 
published a short list of globular clusters with spectroscopy. In addition to globular clusters 
identifications from Kundu & Whitmore (2001), Larsen et al. (2003) globular cluster #1 
agrees with the position of our Src. 59. Sources with candidate globular clusters are noted 
in Table 1. There were no available lists of globular clusters for NGC 4382. 



4.3. X-ray Luminosities and Luminosity Functions 

When converting the source count rates into unabsorbed X-ray (0.3-10 keV) luminosi- 
ties, we assumed that each source was at the distance of the target galaxy. We then used 
the best-fit Chandra X-ray spectrum of the inner (one effective radius for NGC 4365, two 
effective radii for NGC 4382) resolved sources (Tables 3 and 4, row 3). The conversion fac- 
tors were 4.08 x 10 41 ergs cnt^ 1 for NGC 4365 and 3.93 x 10 41 ergs cnt^ 1 for NGC 4382. 
Column 8 of Tables 1 and 2 list the X-ray luminosities in units of 10 37 ergs s^ 1 , and range 
roughly from 1.1 x 10 38 to 2.3 x 10 39 ergs s" 1 for NGC 4365 and 1.1 x 10 38 to 6.2 x 10 39 ergs 
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s- 1 for NGC 4382. 

In Figure 3, we display the cumulative luminosity functions of all resolved sources in 
the S3 field. Each cumulative luminosity function should be the sum of the point source 
(LMXB) population of the target galaxy and the foreground/background population. We fit 
the luminosity function using the same techniques as we have used previously (Sarazin et al. 
2000, 2001; Blanton, Sarazin, & Irwin 2001; Irwin, Sarazin, & Bregman 2002); a single power- 
law, broken power-law, and a cutoff power-law were all used to model the LMXB population. 
The background source population was modeled as discussed in the previous references. For 
NGC 4365, a single power-law fit was acceptable (at ~54% confidence level). However, the 
single power-law fit for NGC 4382 was rejected at the greater than 94% confidence level. 
Attempts at broken power-law fits were more successful [Ax 2 = —6.2 and —19.0, with two 
less degrees of freedom (dof) compared to the single power-law]; however, the power-law 
exponents at high luminosities were steep (~9), so that these models effectively had a cutoff 
above the break luminosity. As a result, we fit the luminosity functions with a cutoff power- 
law model (Ax 2 = —7.3 and —19.3, with one less dof compared to the single power-law): 



where Lx is the X-ray luminosity, L37 = Lx/ (10 37 erg s" 1 ), and L c is the cutoff luminosity. 
The best fits were determined by the maximum likelihood method, and the errors (90% 
confidence interval) were determined by Monte Carlo techniques. For NGC 4365, N Q = 



2.281^ x 10~ 2 , a = lMtoH, and L c = 226lj£ x 10 37 ergs s -1 ; while for NGC 4382, 
N = 0.1912:51, a = 0.93iS:ii, and L c = IIOI23 x 10 37 ergs s -1 . 



At the faint end of the luminosity function, incompleteness may play a role. Since the 
incompleteness is a complicated function of source counts, spectra, background, the point 
spread function, and the detection band, we used simulations to understand the incom- 
pleteness. With MARX 8 simulations, we found that incompleteness can lead to a depressed 
luminosity function for luminosities ~ 30% higher than the minimum detectable luminosity 
in these galaxies. In our cases, that corresponds to ~ 1.4 x 10 38 ergs s -1 (14 counts) for 
both galaxies. This is in rough agreement with other reported simulations (Kim & Fabbiano 
2003). The largest effect of incompleteness for our fits would be in their power- law slopes. 
When we fit the luminosity function at a minimum of 2 x 10 38 ergs s _1 , we find agreement 
within the errors reported above. 




(1) 



8 See http://space.mit.edu/CXC/MARX/. 
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The X-ray luminosity functions of LMXBs have been determined with Chandra in a 
number of other early-type galaxies. In most cases, the observations are limited to Lx > 10 38 
ergs s _1 . The best-fit slope a of NGC 4365 is steeper than the low luminosity slope ol\ in 
NGC 1553 (Blanton et al. 2001), NGC 4649 (Randall et al. 2003), and NGC 4697 (Sarazin 
et al. 2001), although the values for NGC 4365, NGC 4649, and NGC 4697 agree within the 
90% confidence intervals. The slope in NGC 4365 is also flatter than the high luminosity 
slope a h in these three galaxies, and the single power-law slope in NGC 1316 (Kim & 
Fabbiano 2003). The 90% confidence intervals of NGC 4382 a overlap with a { in NGC 1553, 
NGC 4649, and NGC 4697. 

The cutoff luminosities of NGC 4365 and NGC 4382 (« 0.9-3.2 xlO 39 ergs s' 1 ) are 
much higher than the break or cutoff luminosities of NGC 1291, NGC 1553, NGC 4649, and 
NGC 4697 (« 2-6 xlO 38 ergs s" 1 ) (Sarazin et al. 2001; Blanton et al. 2001; Irwin et al. 2002; 
Randall et al. 2003). The differences probably are too large to be explained by distance errors. 
When we fit the luminosity functions in NGC 4365 and NGC 4382 with broken power-laws, 
the break values were somewhat lower than the cutoff luminosities but consistent within the 
errors. Moreover, the break luminosities were significantly larger than the values found in 
NGC 1291, NGC 1553, and NGC 4697. A very high break luminosity or single power-law 
fit was found previously for NGC 1316 (Kim & Fabbiano 2003); this galaxy had a single 
power-law fit, but a break at > 5 x 10 38 ergs s" 1 could not be excluded. 

In Sarazin et al. (2000), we argued that a break luminosity similar to the Eddington 
luminosity of a 1.4M spherically accreting neutron star indicated that the LMXBs above 
the break contain accreting black holes. The large luminosities of the cutoff or broken power- 
laws for NGC 1316, NGC 4365, and NGC 4382 seem less consistent with this interpretation, 
barring a large error in the distance. 

In the luminosity function of NGC 4382, the best-fit clearly misses the brightest source, 
Src. 38. Since Src. 38 is a bright background AGN, this is not indicative of a poor fit. 

4.4. Hardness Ratios 

We determined X-ray hardness ratios for the sources, using the same techniques and 
definitions we used previously (Sarazin et al. 2000, 2001; Blanton et al. 2001; Irwin et al. 
2002). Hardness ratios or X-ray colors are useful for crudely characterizing the spectral 
properties of sources, and can be applied to sources that are too faint for detailed spectral 
analysis. We define two hardness ratios as H21 = (M — S)/(M + S) and H31 = (H — 
S)/(H + S), where S, M, and H are the total counts in the soft (0.3-1 keV), medium (1-2 
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keV), and hard (2-10 keV) bands, respectively. Since the hardness ratios measure observed 
counts, Galactic absorption, as well as QE degradation, must be taken into account when 
interpreting the ratios. The hardness ratios and their la errors are listed in columns 9 and 
10 of Tables 1 and 2. Figure 4 plots H31 vs. H21 for the 48 and 39 sources of NGC 4365 
and NGC 4382, respectively, with at least 20 net counts. The hardness ratio for the sum of 
those sources are (H21,H31) = (0.04, -0.25) for NGC 4365 and (0.02, -0.25) for NGC 4382. 
Sources with ~40 net counts had errors similar to the median of the uncertainties, ~ 0.26. 
The errors scale roughly with the inverse square root of the net counts. 

As was also seen in the bulge of NGC 1291, NGC 1553, and NGC 4697 (Sarazin et al. 
2000, 2001; Blanton et al. 2001; Irwin et al. 2002), most of the sources lie along a broad 
diagonal swath extending roughly from (H21,H31) w (—0.3, —0.7) to (0.4,0.5). In Figure 4, 
the solid line corresponds to Galactic absorption, QE degradation, and power-law spectra 
with photon indices of T = to 3.2. Most of the sources (37/48 and 26/39) are displaced to 
the right of the solid curve. Were the solid curve a good fit to the data, one would expect to 
find roughly equal numbers of sources on each side of the curve. Assuming that this is the 
case, the probability of seeing < 9/48 and 13/39 sources to the left of the line is 0.01% and 
2.66% for NGC 4365 and NGC 4382, respectively. The spectra seem to be more complex 
than a single power-law with varying indices. 

Between NGC 4697 and the bulge of NGC 1291, there were four sources with (H21,H31) 
(—1,-1). These supersoft sources have, essentially, no emission above 1 keV. Based on 
scaling to the 3 sources in NGC 4697, one might have expected ~5-6 supersoft sources 
in NGC 4365 and ~3 supersoft sources in NGC 4382. However, no supersoft sources are 
observed in either galaxy. 

In Figure 4, there are three NGC 4365 sources (Srcs. 36, 69, and 92) and four NGC 4365 
sources with very hard spectra (hardness ratios [H21,H31] > [0.5,0.5]). These may be un- 
related, strongly absorbed AGNs, similar to the sources which produce the hard component 
of the X-ray background, and which appear strongly at the faint fluxes in the deep Chandra 
observations of blank fields (Brandt et al. 2000; Mushotzky et al. 2000; Giacconi et al. 2001). 
There also are two sources with hardnesses ratios of around (—0.4,-1); these are Src. 38 
of NGC 4365 and Src. 10 of NGC 4382; although the hardness ratio of the latter is highly 
uncertain. Studies of other galaxies (Sarazin et al. 2001) and deep blank sky images (e.g., 
Giacconi et al. 2001) suggest that many of sources, similarly lacking hard emission, may 
be unrelated background sources. However, such sources, especially when at smaller radii, 
may be part of the host galaxy. Although the identification of Src. 38 of NGC 4382 with 
a z = 1.19 quasar marks it as a background object, it occupies a typical part of the figure 
([H21,H31] « [-0.2,-0.6]). 



4.5. Variability 



We searched for variability in the X-ray emission of the resolved sources over the duration 
of the Chandra observation using the Kolmogoroff-Smirnov (KS) test (see Sarazin et al. 2001). 
In most cases, the tests were inconclusive. For NGC 4365, six sources (Srcs. 16, 21, 29, 31, 
48, and 86) had a > 95% probability that they were variable. For NGC 4382, only Src. 54 
had a > 95% probability that it was variable. Figure 5 displays the total count histograms 
of the three sources in NGC 4365 (Srcs. 16, 29, and 86) that had a > 99% probability 
that they were variable and NGC 4382 Src. 54 which had a 97% probability that it was 
variable. NGC 4365 Src. 16 seems to be gradually fading. Similarly, NGC 4365 Src. 29 also 
seems to fade gradually, but has only 15 net counts. The X-ray light curve for NGC 4365 
Src. 86 is shown in Figure 6. This source appears to have turned on after the start of the 
observation, and undergone a large outburst after about 1/3 of NGC 4365's exposure. It then 
decreased, and probably underwent at least one more outburst before gradually fading as 
the observation was ending. The light curve of NGC 4382 Src. 54 suggests that it underwent 
an outburst after about 3/8 of NGC 4382's exposure, and that the outburst gradually faded. 
Variable sources are noted in Tables 1 and 2. 



4.6. Spatial Distribution 

As in Sarazin et al. (2001), we have performed two comparisons of the spatial distribution 
of the X-ray sources with that of the optical light; we have fit the position angle (PA) 
distribution and the radial distribution. For both galaxies, we limited the analysis to all 
sources within 3' of the galaxy center to ensure that the circular area lies completely on the 
chip. In this area, we expect ~5 background sources in each exposure. We have adopted 
the Third Reference Catalogue of Bright Galaxies (RC3) values for the optical photometry's 
effective radii, position angle, and ellipticity (de Vaucouleurs et al. 1992). At an effective 
radius of 49"8, NGC 4365 has PA = 40° and e = 0.28; while at an effective radius of 54"6, 
NGC 4382 has e = 0.22. There is no PA listed for NGC 4382 in RC3. Fisher (1997) found 
that NGC 4382 has an isophotal twist, with its PA ranging from 60° at its center to 0° at 
large radii. We adopted that latter PA for NGC 4382. These elliptical isophotes contain 
one-half of the optical light. 

In Figure 7, we show the observed distributions of PAs (modulo 180°, measuring PA from 
north to east). In addition to the observed distribution and the expected optical distribution, 
we have fit a distribution representing the number of expected background sources plus the 
number of sources expected from a set of sources with constant projected density on elliptical 
isophotes. For NGC 4365 and NGC 4382, A% 2 between the optical fits and the maximum 
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likelihood fits indicate that the optical fits are acceptable at the 83% and 54% confidence 
level. The best-fit distribution of NGC 4365 gives PA = 32° and e = 0.29, while the best-fit 
distribution for NGC 4382 has PA = 39° and e = 0.21. The PA value of NGC 4382 lies 
within the range of angles Fisher (1997) reports. 

Adopting the optically determined values for ellipticity and PA, we show the accumu- 
lated source number as a function of semi-major axis, a, in Figure 8. The use of elliptical 
isophotes reduces the number of background sources modeled in each galaxy by ~ 1. We 
modeled the galactic LMXB distribution using a de Vaucouleurs profile, and included a 
uniform distribution of background sources at the level expected from deep Chandra ob- 
servations. If we fix the effective semi-major axis of the LMXB distribution at the value 
determined by the optical light distribution in NGC 4365, a eff = 58"5, the KS test indi- 
cates that the optical distribution is rejected at the > 99% confidence level. The best-fit 
de Vaucouleurs model has an effective semi-major axis of 136". Similarly, the KS test for 
NGC 4382 rejects the optical distribution at the 85% confidence level. The best-fit de Vau- 
couleurs profile for the LMXBs has an effective radius of 129". Both fits indicate that there 
are more X-ray sources at larger radii than expected from the optical light distribution. Ra- 
dial fits from circular isophotes yield similar results. Since NGC 4382 is an SO galaxy, this 
disagreement may reflect the contribution of X-ray sources associated with the disk of the 
galaxy. However, this would not explain the difference between the LMXB distribution and 
the optical light distribution in NGC 4365. 

We note that a significant fraction of the LMXBs in NGC 4365 are associated with 
globular clusters (§ 4.2), as is generally true of elliptical galaxies (Angelini, Loewenstein, 
& Mushotzky 2001; Sarazin et al. 2000). In elliptical galaxies (including NGC 4365), the 
globular cluster population is more broadly distributed than that of the optical light and field 
stars (e.g., Harris 1991). This might help to explain the broader distribution of the LMXBs. 
Figure 8 suggests that there may be two populations of LMXBs, one with a small effective 
radius perhaps consistent with the optical distribution, and one with a larger effective radius. 
It is possible that these two distributions represent the LMXBs formed in globular clusters 
and those formed from field binary star systems. However, we note that in NGC 4697, the 
LMXBs followed the optical light distribution (Sarazin et al. 2001). More elliptical galaxies 
need to be observed to resolve this question. 

To test whether the LMXBs in these galaxies could discriminate between such models, 
we added either a second de Vaucouleurs model or an exponential model to the optical light 
distribution. For either additional component, there were only minimal improvements to 
the fits (Ax 2 > -0.6 for NGC 4365 and > -0.8 for NGC 4382, with one less dof compared 
to the single profile), indicating that the current data cannot statistically discriminate be- 
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tween single de Vaucouleurs distributions with large effective radii or multiple component 
distributions. 



5. Spectral Analysis 

We extracted spectra of the sources and diffuse emission in both galaxies. We restricted 
the spectral analysis to the 0.7-9 keV range. The lower limit was taken to avoid calibration 
uncertainties, while there are few non-background counts beyond 9 keV. We also found a 
possible spectral artifact around 1.6-1.9 keV, which has been noticed by others 9 . We believe 
this artifact is minor enough that it should not effect continuum models (bremsstrahlung 
or power-law); however, we chose to excise the 1.6-1.9 keV band when we fit line models 
(mekal). All of the spectra were grouped to have at least 25 counts per spectral bin prior 
to background correction to enable our use of x 2 statistics. The use of minimum counts per 
spectral bin and restricted energy ranges can result in bins being excluded in the allowed 
energy range. 

The results of the spectral fits are summarized in Tables 3 and 4. Spectra were extracted 
for the resolved point sources ('Sources'), the unresolved diffuse emission excluding the point 
sources ('Unresolved'), and for the total emission (sources and unresolved emission). The 
total emission spectra were studied for comparison to instruments with poorer spatial reso- 
lution that cannot resolve the point sources. The third column gives the geometric region for 
the spectrum; 'Field' implies the entire S3 chip. The value of the absorbing column density 
(Nh) applied to all components of the model emission spectrum is given in column 4. In this 
and other columns, values in parentheses are fixed (not allowed to vary). The fixed value 
of Nh is the Galactic value from Dickey & Lockman (1990). As in Sarazin et al. (2001), 
a two component model was necessary to fit the total spectrum data. A hard component, 
modeled by thermal bremsstrahlung ('bremss') or by a power-law ('power') spectrum, was 
used to model the resolved and unresolved LMXBs. An additional soft line emitting com- 
ponent ('mekal') was used for gaseous diffuse emission. We used the xspec acisabs model 
to correct for the QE degradation. For the hard component, columns 5-7 give the spectral 
model, the temperature Th (for bremsstrahlung) or photon number spectral index T, and 
the unabsorbed flux of the hard component, F\ (0.3-10 keV). Similarly, columns 8-10 give 
the temperature T s , overall heavy element abundance relative to solar, and flux for the soft 
mekal component. For the unresolved emission, the spectra exclude regions around each of 
the resolved sources. The last two columns give the number of net counts in each spectrum, 



See http://asc.harvard.edu/cal/Links/Acis/acis/Cal_projects/index.html. 
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and x 2 per degree of freedom (dof) for the best-fit model. All errors reported in the spectral 
analysis are 90% confidence level errors. Parentheses are used to indicate a frozen parameter 
and brackets are used when an error is unconstrained. 

The background spectra for the resolved sources were determined locally, using the same 
nearby regions as discussed in § 4. For the spectra of the resolved diffuse emission and the 
total spectrum (sources and unresolved emission), the background was taken from the deep 
blank sky backgrounds compiled by Maxim Markevitch 3 . 

The spectra of several spatial regions were analyzed. In NGC 4365, the spectra of 
sources, unresolved emission, and total emission were derived from an inner region, which 
was determined by the elliptical optical isophote containing one-half of the optical light. We 
will refer to this region as within "one effective radius". The semimajor axis of this isophote 
is defined as a e g . In NGC 4365, this isophote has a c g = 58'/5, a semi- minor axis of 42 / /5, and 
a PA of 40° (de Vaucouleurs et al. 1992). In NGC 4382, the corresponding isophote did not 
have sufficient counts from resolved sources for spectral analysis. Therefore, we chose the 
innermost region of NGC 4382 to be twice as large (2 a cS ), which gave an elliptical region 
with a semi-major axis of 123", a semi-minor axis of 95'/5, and a PA of 0° (de Vaucouleurs 
et al. 1992; Fisher 1997). This isophote contains 69% of the optical light. We fit the sources 
(resolved emission), unresolved emission, and total emission in these innermost regions. We 
also examined the spectra of the sources for the entire field. In order to search for changing 
unresolved emission with radius, we chose to analyze the elliptical annulus between one and 
three effective radii for both galaxies, corresponding to ~ 28.8% of the optical light. For 
both galaxies, we used the same PA as the innermost region. In NGC 4365, there were 
enough counts from resolved sources between one and three effective radii to also analyze 
that spectrum. Finally, we analyzed the total emission spectra in the three effective radii 
region. Since the different spectra (source, unresolved, and total) were binned separately, 
all of the counts for a point source were assigned to a region based the source was near the 
edge of the region, and the source spectra used a local background, while the unresolved 
and total spectra used a blank-sky background, the counts of the source spectrum plus the 
unresolved spectrum in a region will not exactly equal the counts of the total spectrum. 

5.1. X-ray Spectra of Resolved Sources 

5.1.1. NGC 4365 

For NGC 4365, we extracted the spectrum of all resolved sources within one effective 
radius. The observed spectrum is shown on the left side of Figure 9. As has been found 
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in other early-type galaxies (Sarazin et al. 2001), the combined spectra of the sources was 
reasonably well-fit by either a thermal bremsstrahlung model with kT^ = 6.08 keV (Table 3, 
row 1) or power-law model with a photon number spectral index of T = 1.67 (row 3). The 
power-law model gave a slightly better fit, so we adopted this model as our best-fit. The 
fits were not improved significantly when the absorbing column was allowed to vary (rows 
2 & 4), so we fixed the hydrogen column at the Galactic value, Nh = 1.63 x 10 20 cm" 2 
(Dickey & Lockman 1990). The adopted best-fit power- law model with Galactic absorption 
and residuals to the fit are shown on the left side of Figure 9. Although the hardness 
ratios indicate that each source is probably fit by a more complex model than a power- 
law and Galactic absorption, the summation of counts from objects with different emission 
properties can be well fit with a single power-law and Galactic absorption. Extending the 
spectrum down to 0.3 keV to see how the restricted energy range used in fitting affects the 
absorbing column does result in a somewhat higher absorption value; however, we believe the 
calibration at low energies is still too unreliable to confidently claim an excess absorption. 

We also searched for any radial variation in the spectra of the sources. We extracted 
the spectrum from an elliptical annulus between one and three effective radii. We first fit 
this region using the source best-fit model from the inner one effective radius region, which 
provided an acceptable fit (row 5). When we allowed the power-law index to vary (row 6), 
we found a better fit (A% 2 = —2.9 with one less dof); however, the two values of T agree 
their 90% confidence limits. Freeing the absorption also yielded a better fit, (A% 2 = — 1.9 
with one fewer dof); however, Galactic absorption was allowed at the 90% confidence limit. 

We also attempted to fit the collective spectrum of all of the resolved sources on the 
S3 chip, which had 2660 net counts. Again, freeing the power-law index from the best-fit 
value for the inner effective radii produced a better fit; however, the power-law index is 
still consistent with the best-fit 90% confidence limit (rows 8 & 9). Freeing the absorbing 
column did not produce a better fit (rows 9 & 10). We find no clear trend between radius 
and spectral fit. This is consistent with results from a study of the sources in 15 galaxies by 
Irwin et al. (2003). 



5.1.2. NGC 4382 

For NGC 4382, the lower number of sources and counts required using the inner two 
effective radii to attain 921 net counts in the spectrum. For this region, the adopted best-fit 
model had a power-law index of 1.52 with the Galactic absorption column of 2.45 x 10 20 cm -2 
(Table 4, row 3). Neither using a bremsstrahlung model (row 1) nor freeing the absorbing 
column (rows 2 and 4) gave a better fit. Again, we remark that although hardness ratios 
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indicate that each source is probably fit by a more complex model than a power-law and 
Galactic absorption, the summation of counts from objects with different emission properties 
can be well fit with a single power-law and Galactic absorption. The best-fit indicated no 
excess absorption when we extended the spectrum down to 0.3 keV. 

Since there were not enough resolved source counts between one and three effective 
radii, we determined the spectrum of the entire field (which of course includes the inner 
two effective radius region discussed above). There were no large statistically-significant 
differences in the fits for the entire field and for the inner region (rows 5-7). Within the 90% 
confidence level, the sources could be fit by a power-law fit with indices between 1.45 and 1.65 
(inclusive of the inner two effective radius fit) and an absorbing column of < 11 x 10 20 cm -2 . 
Note that an important contribution to the spectrum in the outer parts of the observation 
of NGC 4382 is Src. 38, which is a known background AGN. However, this AGN has X-ray 
hardness ratios which indicate that it doesn't have a particularly hard or absorbed spectrum 
(Table 2). The X-ray hardness ratios for this source are similar to the average values for 
the LMXBs in NGC 4382. Since we can rarely identify an AGN with a particular source, 
AGN always contribute to spectra describing the resolved sources. Removing Src. 38 from 
the spectral analysis would hinder comparisons with other galaxies' resolved sources. 

5.2. X-ray Spectra of Unresolved Emission 

5.2.1. NGC 4365 

Again, we began by exploring the spectrum of the inner effective radius, which had 
935 net counts. This spectrum is shown on the left side of Figure 10. First, we attempted 
to model the unresolved emission with a soft mekal component representing the emission 
by diffuse interstellar gas (Table 3, row 11). In general, none of the spectral fits were 
improved by allowing the absorbing column to vary (e.g., row 12), so we will discuss only 
fits with a fixed Galactic column. Values from both mekal-only fits are consistent with 
previous ROSAT measurements from Davis & White (1996). The fits to the unresolved 
emission using only a soft mekal model showed large residuals at high energies. Of course, 
the unresolved emission includes unresolved point sources as well as diffuse gas, and the 
unresolved point sources will have a hard spectrum if they are like the resolved sources. 
Thus, we added a hard component to the fit with the same spectral shape as the best-fit 
model for the source. This produced a dramatic reduction in x 2 by 37 for one less dof, clearly 
indicating values from single temperature fits, like in Davis & White (1996), may be suspect. 
Allowing the power-law index to vary did not improve the fit significantly. Therefore, we 
adopted the model with a power-law index of 1.67, a mekal temperature of 0.561q;o8 keV, and 



-17- 



an abundance in solar units of 0.35 (row 13) as our best-fit model for unresolved emission 
(Figure 10). Although the line emission strength from heavy elements is well constrained, 
minimal emission from hydrogen leads to poorly constrained (> 0.08) abundances. These 
best-fit values are consistent with some of the soft temperatures and abundances found from 
two-temperature fits to ASCA observations (Buote & Fabian 1998). 

We also examined whether the unresolved emission's spectrum changed with radius. The 
unresolved spectrum between one and three effective radii had 1155 net counts. We first fit 
this spectrum with the best-fit model for unresolved emission (row 14). Then, we allowed 
the soft component model to vary in the fit. This produced a much better fit (A% 2 = —37) 
with a statistically significant higher temperature (row 15). This suggests that there is a 
positive radial gradient in the temperature of the gas. If there is a significant gradient, 
this would also affect the spectra in the inner parts of the galaxy. We would be observing 
both cooler gas located near the center of the galaxy, and hotter gas at larger radii seen in 
projection against the center. Thus, we tried fitting the emission from the inner effective 
radius with two soft components and a hard component; that model did not improve on the 
single soft+hard model spectral fit we report in Table 3. A temperature gradient may exist 
in the diffuse gas; however, deeper observations of NGC 4365 are necessary to do a more 
detailed deprojection. 

In addition to a possible temperature gradient in the diffuse emission, the spectral fits 
of the sources and unresolved emission point to a possible change in the unresolved sources 
with respect to radius. In the resolved sources, the model flux for sources in the inner 
effective radius and sources between one and three effective radii are about the same. The 
correction for increased contribution by background sources based on the luminosity function 
we fit earlier suggests that the flux from resolved sources between one and three effective 
radii should be ~ 0.9 times the flux from resolved sources in the inner effective radius. The 
luminosity function of background sources in a particular field can differ from an average field 
either by total number of sources or by having a brighter than expected background source 
in the field. Source 69 is a very hard and bright source that might be such an unusually 
bright AGN. This source accounts for ~ 1/4 of the hard counts in the one to three effective 
radii resolved source region. If we were to exclude this source, the flux ratio between the two 
regions would be ~ 0.65. One would expect that if the unresolved sources distribution is the 
same as the resolved source distribution, we should find a flux ratio between ~ 0.6 and 1.0 
for the unresolved sources. Instead, the best-fit unresolved emission between one and three 
effective radii is a purely soft component fit. At best, we can limit the ratio to be < 0.15. 
This suggests that the unresolved sources radial distribution may decline more rapidly than 
the resolved source distribution. Poisson errors in the 2-9 keV count rate of the one to three 
effective radii unresolved emission are too large to statistically confirm this result. Deeper 
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X-ray observations of NGC 4365 are needed to resolve enough sources to allow the faint 
source and bright source radial distribution to be compared. 

5.2.2. NGC 4382 

The spectrum of the unresolved emission of the inner two effective radii (right side of 
Figure 10) had 2620 counts. We initially fit the spectrum with a model containing only a soft 
mekal component representing emission from diffuse interstellar gas (row 8). The best-fit 
models with varying absorption gave values below the Galactic column, but with an upper 
limit which included the Galactic value (Table 4, row 9). Thus, we chose to use the Galactic 
column for all remaining fits of this galaxy's spectra. As with NGC 4365, the addition of 
a hard component representing unresolved sources led to a much better statistical fit (row 
10). Allowing the parameters of the hard component to vary did not result in a significantly 
better fit, and the hard component parameters were very poorly constrained. Thus, we fixed 
the shape of the spectrum of the hard component to the adopted model for the resolved 
sources in the galaxy. This leads to an acceptable fit (row 10); however, its abundance is 
abnormally high and poorly constrained due to minimal emission from hydrogen. Therefore, 
we have frozen the abundance at 0.5 solar for all remaining soft component fits in this galaxy. 
This choice (row 11) is very close to the best fit (Ax 2 = 0.3). This best-fit model has a 
temperature of 0.28^03 

and is shown in Figure 10. 

We extracted the unresolved emission spectrum from a larger (but overlapping) elliptical 
annulus extending from one to three effective radii. The spectrum in this region was reason- 
ably fit by the same gas temperature and abundance which provided the best fit within two 
effective radii (row 12). When the gas temperature was allowed to vary (row 13), the best-fit 
temperature increased; however, that increase was not statistically significant. Unlike in 
NGC 4365, we find no change in the flux ratio of the resolved and unresolved sources with 
increasing radius. 

5.3. Total X-ray Spectra 

5.3.1. NGC 4365 

We first extracted the spectrum of all of the emission within one effective radius. The 
spectrum in this region is consistent with the best-fit spectra for the resolved and unresolved 
emission (row 16); when the spectral parameters of the hard and soft components were al- 
lowed to vary, the fits were not significantly improved (rows 17 and 18) and did not constrain 
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the abundance. We also attempted a model with a bremsstrahlung hard component (row 
19). This fit slightly better than the power-law model, was more successful at constrain- 
ing the abundance (< 0.26 solar), but did not constrain the bremsstrahlung temperature 
well (> 5.21 keV). We therefore adopted the sum of the best-fit spectra for the resolved 
and unresolved emission as our best-fit for the total emission. The total emission within 
three effective radii yielded the same fit within the overlap of the two region's errors for a 
power-law hard component model. Again, a bremsstrahlung hard component could better 
constrain the soft component abundance at the expense of constraining the hard component 
temperature. The abundance is significantly lower than the best-fit unresolved abundance; 
however, the statistics can not discriminate between the power-law fit, with its poorly con- 
strained abundance (> 0.01 solar), and the bremsstrahlung fit to within a 90% confidence 
level. Therefore, we do not claim any abundance gradient in the total emission. 

In Matsumoto et al. (1997), a bremsstrahlung model with a temperature of 12. O^ 2 ^ 
keV or a power-law model with index 1.8^q;4 was fit to the hard component, while the soft 
component was fit by either a mekal or a Raymond-Smith model with temperature ~ 0.5 
keV and poorly constrained abundance (> 0.016 solar). These results are consistent with 
what we have found. 



5.3.2. NGC 4382 

The results for NGC 4382 were similar to those for NGC 4365. Within two effective 
radii, the total spectrum was fit well by the combination of the best-fit hard component 
from the resolved sources and the best-fit soft component from the diffuse emission (Table 4, 
row 14). Freeing the hard and soft spectral parameters did not improve the fit significantly 
(rows 15 and 16). A bremsstrahlung hard component did not produce a better fit and poorly 
constrained the bremsstrahlung temperature > 10.38 keV. There was no significant change 
between the total emission within two effective radii and the total emission within three 
effective radii. 

Once again, the results were similar to Matsumoto et al. (1997). They used the same 
hard component model discussed above for NGC 4365, and fit a soft temperature ~ 0.3 keV. 
Matsumoto et al. (1997) were able to constrain the abundance using the Raymond-Smith 
model to 0.029-0.21 solar, but could not constrain the abundance for the mekal model. 
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5.4. Discrete versus Diffuse Emission 

We used both counts (0.3-10.0 keV) and spectral fitting results to determine the fraction 
of the X-ray emission which was due to diffuse gas versus the fraction due to point sources, 
whether resolved or not. The unresolved emission is a combination of unresolved point 
sources and truly diffuse gaseous emission. We used the resolved source hardness ratios of 
the each region to separate unresolved counts into point sources and diffuse gas. We also 
compared the hard flux (point sources) to the soft flux (diffuse gas) in the total emission 
spectra. Uncertainty in the spectral models used lead to a much wider range of derived 
flux ratios than count ratios. Since the gas is softer than the sources, absorption affects the 
diffuse emission more. Therefore, the unabsorbed flux ratio of point sources to total emission 
need not be larger than the count ratio of point sources to total emission, as one one would 
expect. 

For the inner one effective radius of NGC 4365 the point sources are responsible for 
75 ± 7% of the total counts and ~ 70% of the unabsorbed flux (Table 3, rows 16-19). 
Approximately 60% of the point source counts are resolved. The unresolved point source 
flux overlaps with 90% confidence level fluxes obtained from the derived luminosity function 
extended down to 10 36 ergs s _1 . Moving to the region within three effective radii the point 
sources are responsible for 75 ± 6% of the total counts and ~ 55% of the unabsorbed flux 
(rows 20 & 21). 

In the inner two effective radii of NGC 4382 the point sources account for 55 ± 5% of 
the total counts and ~ 50% of the unabsorbed flux (Table 4, rows 14-17), with ~ 40% of 
the point source counts resolved. Again, the unresolved point source flux overlaps with 90% 
confidence level fluxes obtained from the derived luminosity function extended down to 10 36 
ergs s _1 . In the inner three effective radii, the point sources account for 56 ±6% of the total 
counts and ~ 50% of the unabsorbed flux (rows 18 & 19). 

For both galaxies, it was crucial that a majority of the source were resolved. For power- 
law models, the 90% confidence limit on the power-law index was cut by a factor of 2-3 
sources were resolved. If one uses bremsstrahlung models, the bremsstrahlung temperature 
had no upper constraint unless the sources' spectra could be fit separately. 

6. Spatial Distribution of the Diffuse Gaseous Emission 

We derived the radial distribution of the diffuse gas, based on the surface brightness 
distribution in the soft band (0.3-1 keV). We excluded regions around each of the resolved 
sources in determining the diffuse surface brightness. Although the diffuse gas generally 
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dominates in the soft band, this emission still contains a contribution from unresolved sources 
(§ 5.2). To remove the effect of the unresolved point sources, we have assumed that the 
spectrum of the sum of all resolved sources in the S3 field is representative of the spectrum 
of the unresolved sources. Based on the observed spectra, we also assume that the diffuse 
gas produces no significant hard band (2-10 keV) emission. Thus, we determined the surface 
brightness profile in the hard band (removing resolved sources), and scale this by the ratio 
of soft to hard band emission expected based on the best-fit spectrum of the resolved point 
sources. This soft band surface brightness profile for the sources was subtracted from the 
observed soft band profile, and we attribute the remainder to diffuse gas. We performed this 
operation in regular circular annuli out to ~ 3', and corrected the profiles for background 
and exposure. 

Figure 11 displays the observed soft band gaseous surface brightness profiles as a function 
of projected radius, r. The dashed lines display the best-fit de Vaucouleurs profiles with 
effective radii fixed at the values for the optical light in the galaxies (49 / /8 for NGC 4365 
and 54'/ 6 for NGC 4382). The normalization was varied to achieve the best fit. NGC 4365 is 
acceptably fit by such a profile, with a x 2 of 29.8 for 29 dof. NGC 4382 is unacceptably fit by 
the de Vaucouleurs profile, and the diffuse gas surface brightness is consistently larger than 
the optical profile at r > 16". Baggett, Baggett, & Anderson (1998) found that NGC 4382 
could be decomposed into a bulge with an effective radius of 71.4" and a disk with a scale 
length of 165.5". Adjusting the normalization of the bulge plus disk profile better fit the 
X-ray profile; however, it was still an unacceptable fit (x 2 > 99 for 29 dof). Therefore, we 
attempted a standard beta model profile, 



where r c is the core radius. The fit for NGC 4365 was much better than the de Vaucouleurs 
profile, x 2 of 17.7 for 27 dof with r c = (0.50 kpc) and (3 = 0.3981°;°^ (90% confidence 

errors). This near power-law fit is similar to those for NGC 1291 and NGC 4697 (Irwin et al. 
2002; Sarazin et al. 2001). NGC 4382 had a similar quality beta model fit with r c = 35"5± 2 3^ 
(3.2 kpc) and (3 = 0.599 ± 0.265. Since /3 is constrained by the outer data points, the much 
larger core radius of NGC 4382 is responsible for the larger confidence limit of its (3 compared 
to NGC 4365. The beta models for both galaxies are similar to the profiles of X-ray bright 
galaxies within the errors (Forman et al. 1985; Trinchieri et al. 1986). For NGC 4382, its 
large core radius, rounder X-ray emission than optical emission, and lower gas temperature 
at inner radii are roughly consistent with predictions by Brighenti & Mathews (1996) for a 
rotating interstellar gas distribution; however, the core radius is about half their predicted 
value for a galaxy in Virgo. 
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7. Conclusions 

The ability of Chandra observations to separate X-ray emission into point sources 
(LMXBs and background sources) and diffuse emission (unresolved LMXBs and diffuse 
gas) plays a crucial part in exploring the properties of the X-ray faint early-type galax- 
ies NGC 4365 and NGC 4382. We detected 99 sources in NGC 4365 and 58 sources in 
NGC 4382, while we expect that approximately f 2 sources in each galaxy are due to unre- 
lated background sources. Within one effective radius of NGC 4365, ~ 45% of the counts is 
resolved into sources, ~ 30% of the counts is attributed to unresolved LMXBs, and ~ 25% 
is attributed to diffuse gas. This distribution is consistent with that of another X-ray faint 
elliptical galaxy, NGC 4697 (Sarazin et al. 2001), after correcting for different sensitivities. 
Within two effective radii of NGC 4382, ~ 22% of the counts is resolved into sources, ~ 33% 
of the flux is attributed to unresolved LMXBs, and ~ 45% is attributed to diffuse gas. The 
higher diffuse gas fraction in this lenticular is consistent with NGC 1553 (Blanton et al. 
2001), also an X-ray faint lenticular galaxy; however, NGC 4382 resolves less of its flux into 
sources than NGC 1553, despite a better sensitivity. 

The hardness ratios of the resolved sources tend to be harder (more H21) than a range of 
power-law indices with Galactic absorption. This suggests that the source emission spectra 
are more complex than single power-laws. No supersoft sources were found in either of 
the two galaxies. The spectra of the sum of all the sources were best fit by a power-law 
model with Galactic absorption. The power-law of NGC 4382, 1.52, was harder than that 
of NGC 4365, 1.67; however, the indices are the same within their 90% confidence interval 
and are consistent with the best-fit value of the sources simultaneously fit in a survey of 15 
galaxies (Irwin et al. 2003). 

For sources in NGC 4365, 18 sources out of 37 in a WFPC2 field of view had positions 
correlated with a globular cluster compared to the ~ 0.6 expected from random association. 
An additional source's position has also been identified with a globular cluster. Three of the 
sources in NGC 4365 are variable at the >99% confidence level. The spatial distributions of 
the sources are broader than the optical de Vaucouleurs distribution of each galaxy. 

The most surprising results from analysis of the sources are the luminosity fits. Both 
galaxies were best fit with cutoff power-laws, where the cutoff luminosities were between 
~ 0.9-3.1 xlO 39 ergs s _1 . Although broken power-laws were within the 90% confidence 
interval, the break luminosity was similar to the cutoff luminosity. This is in contrast to 
the prior break/cutoff luminosities of NGC 1291, NGC 1553, and NGC 4697 (« 2-6 xlO 38 
ergs s _1 ), which had been identified with the Eddington luminosity of a 1.4 M spherically 
accreting neutron star (Sarazin et al. 2001; Blanton et al. 2001; Irwin et al. 2002). It seems 
unlikely that the cutoff luminosities of NGC 4365 and NGC 4382 can be attributed to that 
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physical mechanism. One other X-ray faint early-type galaxy, NGC 1316, indicates a break 
luminosity that may be as high as the cutoff luminosities found in NGC 4365 & NGC 4382; 
however, the galaxies in this paper provide a tighter constraint on a cutoff/break luminosity 
well above the neutron star Eddington luminosity. 

In NGC 4365, the temperature of the diffuse gas is ~ 0.6 keV and its abundance is 
poorly constrained (> 0.08 solar) in the inner effective radius. There is some evidence for a 
positive temperature gradient in NGC 4365; however, deeper observations are necessary to 
confirm any radial gradients in temperature or abundance. NGC 4382 was best fit with a 
temperature of ~ 0.28 keV and a poorly constrained abundance (> 0.13 solar). 

The diffuse gas surface brightness profiles in both galaxies were fit by standard beta 
model profiles. NGC 4365 had a small (5", 0.5 kpc) core radius and f3 ~ 0.40, similar to 
NGC 4697. NGC 4382 had a larger core radius ~ 36" (3.2 kpc) and (3 between ~ 0.33 
and 0.86. The larger core radius, combined with the X-ray emission being rounder than the 
optical emission, and evidence of a cooler central region, may indicate rotation of the diffuse 
gas. 

We are very grateful to Arunav Kundu for providing us with his unpublished list of 
globular clusters in NGC 4365, and for several very helpful conversations. We thank Mike 
Skrutskie for his help and advice in using the 2MASS catalog. Support for this work was 
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Numbers GO1-2078X and GO1-3100X, issued by the Chandra X-ray Observatory Center, 
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under contract NAS8-39073. This research has made use of the NASA/IPAC Extragalactic 
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of Technology, under contract with the National Aeronautics and Space Administration, the 
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produced at the Space Telescope Science Institute under U.S. Government grant NAG W- 
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Table 1: Discrete X-ray Sources in NGC 4365 
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85 


7 


20 


28 


2 


99 


9 


2 


81+0 


88 


3 


20 


11 


50 


-0 


25±S 


49 
39 


-1 


°°±o 


00 
00 




58 


CXOU 


J122425 


1+071734 


12 


24 


25 


13 


7 


17 


34 


9 


100 


4 


10 


92+1 


70 


6 


42 


44 


58 


+0 


4112 


20 
24 


+0 


ll+° 


29 
31 




59 


CXOU 


J122424 


0+071745 


12 


24 


24 


04 


7 


17 


45 


1 


100 


7 


3 


35+0 


96 


3 


49 


13 


67 


-0 


33^2 


57 
40 


-0 


"IS 


47 

42 


d,f 


60 


CXOU 


J122431 


1+071724 


12 


24 


31 


18 


7 


17 


24 


7 


108 


2 


7 


62+1 


42 


5 


35 


31 


11 


+0 


02+° 


30 
30 


+0 


0312 


31 
31 
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TABLE 1 — Continued 



Src. 
(1) 




Name 
(2) 




R.A. 

(h:m:s 
(3) 


) 


Dec. 

(":':": 
(4) 




d 
(") 
(5) 


Count Rate 
CIO" 4 s- 1 ) 
(6) 


SNR 

(7) 


(8) 




(9) 






(10) 


61 


CXOU 


.1122424 


9+072046 


12 


24 


24 


92 


7 


20 


46 


9 


114 


.4 


2 


,87±0. 


95 


3 


03 


11 


73 


-0 


,08+g 


52 


-0 


.27+0 


62 


CXOU 


J122435. 


9+071925 


12 


24 


35 


.93 


7 


19 


25 


.9 


116 


.3 


3 


.82+1. 


03 


3. 


72 


15 


.59 


+0 


00+u 


'II 

.83 


+0 


69+0 

— U 


63 


CXOU 


J122420. 


7+071938 


12 


24 


20 


.72 


7 


19 


38 


.0 


117 


.2 


6 


.59+1 


34 


4. 


93 


26 


.93 


+0 


.12+° 


30 
.33 


-0 


ll+o 


64 


CXOU 


J122421. 


6+071757 


12 


24 


21 


.62 


7 


17 


57 


.6 


118 


.8 


3 


.94+1. 


04 


3. 


80 


16 


.10 


+1 


00+u 


00 


+1 


00+0; 


65 


CXOU 


J122435. 


3+071805 


12 


24 


35 


.30 


7 


18 


05 


.4 


119 


.9 


10 


.81 + 1 


70 


6. 


34 


44 


13 


-0 


08+° 


'21 
.20 


-0 


.71+0 


66 


CXOU 


J122425. 


8+071706 


12 


24 


25 


.81 


7 


17 


06 


.1 


123 


.2 


3 


.41+0. 


96 


3. 


54 


13 


.92 


-0 


33+0 


46 


-0 


66+;; 


67 


CXOU 


J122423. 


3+071720 


12 


24 


23 


.35 


7 


17 


20 


.6 


126 


.4 


3 


.60+0 


99 


3. 


62 


14 


.70 


+0 


29+0 


42 


+0 


.17+0, 


68 


CXOU 


J122435. 


8+071806 


12 


24 


35 


.88 


7 


18 


06 


.2 


127 


.2 


3 


.87+1. 


04 


3. 


73 


15 


.82 


+0 


31+0 


.38 


+0 


08+0; 


69 


CXOU 


J122426. 


3+071653 


12 


24 


26 


.39 


7 


16 


53 


.2 


133 


.5 


48 


.43+3 


54 


13. 


67 


197 


.75 


+0 


65+0 


'08 


+0 


64+0- 


70 


CXOU 


J122426. 


7+071650 


12 


24 


26 


.70 


7 


16 


50 


.2 


135 


.6 


8 


.21 + 1. 


47 


5. 


58 


33 


.52 


+0 


24+0 


24 


-0 


.16+0 


71 


CXOU 


J122420. 


1+071739 


12 


24 


20 


.18 


7 


17 


39 


.4 


146 


.8 


16 


.46+2 


08 


7. 


89 


67 


.20 


+0 


04+0 


.16 


-0 


.44+0 


72 


CXOU 


J122437. 


0+071757 


12 


24 


37 


.08 


7 


17 


57 


.3 


147 


.1 


3 


.59+1 


01 


3. 


57 


14 


.65 


+0 


46+0 


30 


-0 


.29+0 


73 


CXOU 


J122423. 


9+071651 


12 


24 


23 


.95 


7 


16 


51 


.0 


147 


.4 


5 


.06+1 


17 


4. 


34 


20 


66 


+0 


.47+0 


'26 


-0 


.33+0 


74 


CXOU 


J122438. 


4+071902 


12 


24 


38 


.46 


7 


19 


02 


.8 


151 


.8 


4 


.67+1. 


21 


3. 


86 


19 


.05 


-0 


.10+0 


49 


+0 


.15+0, 


75 


CXOU 


.1122419 


2+072019 


12 


24 


19 


.25 


7 


20 


19 


.3 


153 


.9 


2 


.88+0. 


90 


3. 


20 


11 


78 


+0 


38+0 


50 


+0 


55+g 


76 


CXOU 


J122418. 


0+071839 


12 


24 


18 


.01 


7 


18 


39 


.1 


154 


.4 


3 


.36+0 


97 


3. 


46 


13 


.70 


-0 


28+0 


'47 


-1 


.00+ 2 „ 

— 


77 


CXOU 


J122418 


0+071803 


12 


24 


18 


.04 


7 


18 


03 


g 


163 


.6 


3 


.45+0. 


98 


3 


52 


14 


08 


+0 


oo±° 


'45 


-0 


ll+o 


78 


CXOU 


J122422 


1+071646 


12 


24 


22 


.19 


7 


16 


46 


1 


164 


.7 


9 


.95+1 


62 


6 


14 


40 


63 


+0 


15+0 


24 


+0 


io+°, 


79 


CXOU 


J122434 


9+072116 


12 


24 


34 


.93 


7 


21 


16 


1 


165 


.4 


2 


.86+0 


95 


3. 


02 


11 


.67 


+0 


09+0 


55 
61 


+0 


.17+°' 


80 


CXOU 


J122423 


4+072136 


12 


24 


23 


.44 


7 


21 


36 


.5 


168 


1 


2 


.77+0. 


91 


3. 


03 


11 


.29 


+0 


.07+0 


47 

.50 


-0 


.37+0 


81 


CXOU 


J122436 


4+071706 


12 


24 


36 


.47 


7 


17 


06 


.7 


169 


.2 


3 


.34+0 


96 


3. 


49 


13 


63 


-0 


.27+0 


51 


-0 


.56+0 


82 


CXOU 


J122417. 


9+071741 


12 


24 


17 


.91 


7 


17 


41 


.4 


174 


.7 


4 


.85+1 


15 


4. 


23 


19 


.78 


+0 


38+0 


.37 


+0 


51+0, 


83 


CXOU 


J122417. 


8+071715 


12 


24 


17 


.80 


7 


17 


15 


.6 


189 


.5 


7 


.11 + 1. 


38 


5. 


17 


29 


.02 


+0 


40+0 
— u 


'25 
.33 


+0 


03+0- 


84 


CXOU 


J122436 


6+072129 


12 


24 


36 


.65 


7 


21 


29 


.3 


191 


.8 


9 


.43+1. 


81 


S. 


22 


38 


.51 


-0 


■n±S 


29 
27 


-0 




85 


CXOU 


J122425 


1+072211 


12 


24 


25 


.19 


7 


22 


11 


.2 


192 


.9 


3 


.33+1 


03 


3. 


24 


13 


.59 


+0 


39t°o 


34 
.49 


+0 


18l° 


86 


CXOU 


J122438 


3+072107 


12 


24 


38 


.31 


7 


21 


07 


.0 


193 


.7 


28 


.87+2 


88 


10 


01 


117 


.88 


-0 




12 
.11 


-0 


«±°0 


87 


CXOU 


J122427 


2+072231 


12 


24 


27 


.28 


7 


22 


31 


.2 


207 


.9 


8 


.82+1. 


87 


4 


72 


36 


.03 


+0 


■32i° 


34 
.44 


+0 


.321?, 


88 


CXOU 


J122414 


7+072004 


12 


24 


14 


74 


7 


20 


04 


.9 


210 


.3 


37 


.47+3. 


20 


11. 


72 


153 


.01 


-0 


22+0 


.11 

.10 


-0 




89 


CXOU 


J122424 


8+072229 


12 


24 


24 


.83 


7 


22 


29 


.4 


211 


.9 


4 


.53+1 


35 


3 


36 


18 


.48 


+0 


73+; 


24 
.10 


+0 


76+;- 


90 


CXOU 


J122441 


0+072047 


12 


24 


41 


05 


7 


20 


47 


.8 


216 


.9 


4 


.95+1 


38 


3 


59 


20 


21 


+0 


86+; 


14 

.67 


+0 


85+; 


91 


CXOU 


J122438 


9+071625 


12 


24 


38 


96 


7 


16 


25 


.8 


224 


.3 


12 


.21 + 1. 


84 


6. 


65 


49 


.87 


+0 


05+0 
uo -0 


20 
.20 


-0 




92 


CXOU 


J122438 


9+072152 


12 


24 


38 


.93 


7 


21 


52 


.4 


231 


.6 


7 


.80+1 


61 


4. 


83 


31 


.85 


+0 


•«±° 


28 
.52 


+0 


68i° 


93 


CXOU 


J122425 


4+072258 


12 


24 


25 


.42 


7 


22 


58 


.7 


238 


.6 


7 


.47+1 


62 


4. 


60 


30 


.48 


+0 


.30±° 


26 
.32 


-0 


-33t° 


94 


CXOU 


J122435 


9+071533 


12 


24 


35 


94 


7 


15 


33 


.3 


239 


.6 


5 


.70+1. 


30 


4. 


39 


23 


.26 


+0 


■n±S 


34 
.37 


-0 


.05+.°, 


95 


CXOU 


J122420 


8+072237 


12 


24 


20 


.84 


7 


22 


37 


.0 


240 


1 


6 


.85+1. 


55 


4. 


41 


27 


96 


+0 


-sot , 


25 
.38 


-0 


ioi°o 


96 


CXOU 


J122440 


9+072204 


12 


24 


40 


.90 


7 


22 


04 


.1 


260 


.6 


55 


.42+3 


99 


13. 


88 


226 


.30 


-0 


19+0 


08 
.08 


-0 


•«±S 


97 


CXOU 


J122437 


9+072311 


12 


24 


37 


.91 


7 


23 


11 


.0 


285 


.9 


22 


.69+2 


74 


8. 


27 


92 


.64 


-0 


■oil? 


16 

.16 


-0 


25+0 
^°-0 


98 


CXOU 


J122419 


9+072332 


12 


24 


19 


91 


7 


23 


32 


.6 


296 


.1 


3 


.83+1 


27 


3 


02 


15 


64 


+ 




53 
.58 


+0 




99 


CXOU 


J122416 


0+072324 


12 


24 


16 


03 


7 


23 


24 


.3 


317 


.7 


36 


.75+4 


33 


8 


65 


150 


05 


+0 


49tg 


.15 
18 


+0 


431° 



Notes 
(11) 



.00 
.37 
.39 
.18 
.74 
.27 
.49 
.59 
.51 
.56 



.44 
.36 
.82 



.46 
.56 
.15 



.08 
.15 
.67 



.55 
.10 
.09 



Note. — The units for Lx are 10 37 ergs s 1 in the 0.3-10 keV band. 

a Sources near the center may be confused with nearby sources, making their positions, fluxes, and extents uncertain. 
''Source is noticeably more extended than PSF. 

c Source is at the edge of the S3 detector, and flux is uncertain due to large exposure correction. 

^Possible faint optical counterpart. 

e Idcntificd with known AGN. 

■^Globular cluster is possible optical counterpart. 

9 Source may be variable. 
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Table 2: Discrete X-ray Sources in NGC 4382 (M85) 



Src. 










R.A. 






Dec. 




d 




Count Rate 




















No. 




Name 




| 


h:i 




) 


( 




:") 




(" 


) 


(10 


-4 - 


) 


SNR 


L - 


K 




H21 






H31 


(1) 




(2) 






(3) 






(4) 




(5) 












7 ) 


(8) 




(9) 






(10) 


1 


cxou 




n i i Qi i 07 


12 


2. r ) 


24 


06 


18 


11 


27 


.7 


1 


.8 


22 


79±2 


55 


8 


93 


89 


.66 


— 


16_ 


13 

.13 





, fi +u. 

36 -0 


2 


cxou 




A 1 1011O1 

. U+lnllZl 


12 


25 


23 


99 


18 


11 


21 


7 


4 


. 2 


16 


44±2 


.15 


7 


.63 


64 


.71 


— 


1 n+° 


17 
.16 


— 0. 


OQ + 0. 


3 


cxou 


,J LZZbZo 


A I 101101 

.9+1811..11 


12 


25 


23 


90 


18 


11 


31 


1 


5 


.5 


10 


90±1 


.79 


6 


09 


42 


.88 


— 


08_ 


25 
24 


— 0. 


68_ 


4 


cxou 




r; i ion ojs 
.iJ-(- 1B1 IZD 


12 


25 


23 


.58 


18 


11 


26 


.7 


6 


(i 


11 


87±1 


88 


6 


33 


46 


.71 


— 


04+°- 
04 -0 


27 

.27 


— 0. 


n7 +0. 
07 -0 


5 


cxou 




C 1 1 O 1 1 O A 

-0 + 181 1..-)4 


12 


25 


24 


54 


18 


11 


34 


.9 


11. 


5 


11 


97±1 


.84 


6 


.52 


47 


.11 


+0 


.35_„ 


24 
.29 


+0. 


oo + O. 
28 -0. 


6 


cxou 


,J izzbZo 


.u-\- lol 1..1D 


12 


25 


23 


(id 


18 


11 


36 


4 


12 


2 


5 


10±1 


.24 


4 


.12 


20 


.06 


+0 


' 26 -0 


37 
.46 


+0. 


1S + 0. 
18 -0. 


7 


cxou 


,J LZZbZo 


.7+ 181 139 


12 


25 


23 


. 75 


18 


11 


39 


.9 


14. 


6 


9 


.64±1 




5 


.79 


37 


.92 


— 


20_ 


20 
.18 


— 0. 


65 -0 


8 


cxou 




■"7 1 1 Q1 1 ^1 

, ( -|- lo 1 141 




2. r ) 


24 


70 




11 


41 


(1 




7 


11 


31±1 


.79 


6 


.33 


44 


.48 


+0 


.69_„ 


17 

.31 


+0. 


46 -0. 


9 


cxou 


T 1 OOKOK 
,I IZZbZb 


^ I 1 Ol 1 /|1 

, o+ lo 1 141 


12 


25 


25 


. (i 5 




11 


41 


.5 


2- 




5 


,24±1 


.25 


4 


.20 


20 


.63 


— 


09_ 


45 
.41 


+ 0. 


99+0. 
22 -0 


1(1 


cxou 


7 1 OOr,00 

,I IZZbZZ 




12 


25 


22 


.84 


18 


10 


53 


6 


36. 


' 


7 


19±1 


43 


5 


.01 


28 


.28 


— 


A1 -0 


27 
.21 


— 0. 


07+ 1 
97 -0 


1 1 


cxou 




.8+181(148 


12 


25 


24 


.83 


18 


10 


48 


,7 


38. 


9 


3 


.32±0 


.99 


3 


.36 


13 


.08 


— 




63 
.41 


— 0. 


21 -0 


12 


cxou 


J122524 


.9+181205 


12 


25 


24 


.95 


18 


12 


05 


.7 


41. 


8 


27 


91±2 


.76 


10 


12 


109 


.83 


+0 


1 1 +0- 


14 
.15 


+0. 


nn+0. 
00_ Q 


13 


cxou 


J122526 


.9+181149 


12 


25 


2(1 


94 


18 


11 


49 


,7 


47. 


7 


24 


.71±2 


.59 


9 


.55 


97 


.24 


— 


03„ 


11 
14 


— 0. 


45 -0 


14 


cxou 


J IZZbZu 


. o + 181 Izl 


12 


25 


20 


.53 


18 


11 


21 


.8 


50. 


1 


2 


.74±0 


.91 


3 


.01 


10 


.78 


— 


A1 -0 


86 
47 


+0. 


14 + °. 


15 


cxou 


J IZZbZZ 


1 i i o i a a o 
. 1 + 181(14..-) 


12 


25 


22 


.15 


18 


10 


43 


.5 


50 


3 


17 


38±2 


.22 


7 


.82 


68 


.40 


+0 


01+ 

21_ 


IS 
.19 


— 0. 


O9 + 0- 


16 


cxou 


t 1 oo r on 

J IZZbZu 


,6+iolUbz; 


12 


25 


20 


69 


18 


10 


52 


.5 


58. 


3 


2 


.80±0 


.93 


3 


.01 


11 


.03 


— 


n+0 

n -o 


68 
59 


+0. 


17+0. 


17 


cxou 


7 1 OO C OT 

,J IZZbZ 1 


.4+ lSl lo f 


12 


25 


27 


.47 


18 


11 


57 


.8 


58. 


4 


10 


08±1 


.67 


6 


.03 


39 


.65 


+0 


22_ 


2d 
27 


+0. 


n«+°. 
06 -0 


IS 


cxou 


t 1 ook 1 n 


T 1 1 SI 1 OK 

. / -(- lOl Iz.j 


12 


25 


19 


. 75 


18 


11 


25 


(1 


61. 


2 


4 


.49±1 


.14 


3 


.92 


17 


.66 


— 


07 -0 


38 
.36 


— 0. 


co + 0. 
53 -0 


19 


cxou 


t 1 oo r ot 
,J IZZbZ I 


O 1 1 OIOI c 

.2 + 181zlb 


12 


25 


27 


.29 


18 


12 


15 


.4 


67. 


8 


7 


,24±1 


.56 


4 


.65 


28 


.48 


— 


1/1+° 


32 
.30 


— 0. 


,,+0. 
31_ 


2(1 


cxou 


7 1 OOKOI 

,J lZZoZl 


. ( -\- lolUZO 


12 


25 


21 


. 72 


18 


10 


20 


.3 


68. 


1 


9 


07±1 


.57 


5 


.79 


35 


.69 


+0 


09_ 


28 
.29 


+0. 


na+°. 
09 -0 


21 


cxou 


7 1 oo c o o 
J IZZbZa 


O 1 1 O 1 O A E 

.3+ iolzUo 


12 


25 


28 


.35 


18 


12 


05 


7 


73. 


2 


6 


,11±1 


40 


4 


35 


24 


.04 


+ 1 


oo_ 


00 
.29 


+ !■ 


nn+0. 


22 


cxou 


,J LZZoZl 


1 1 SI OQO 

, D+ lo 1 


12 


25 


21 


. 68 


18 


12 


32 


1 


74. 


2 


3 


.54±1 


.01 


3 


.51 


13 


.91 


— 


,33_ 


57 
.40 


+0. 


n-,+°- 


23 


cxou 


t 1 oo r 1 o 

,J IZZblo 


.4+ lol 1,3 / 


12 


25 


18 


42 


18 


11 


37 


(1 


80. 


9 


14 


54±1 


.98 


7 


.36 


57 


.22 


— 


on+° 
20_ 


22 
.20 


— 0. 


ii+0 

n -o 


24 


cxou 


T 1 OO C 1 A 


O 1 1 O 1 o o c 

,8+lolz3b 


12 


25 


19 


.82 


18 


12 


35 


.3 


91. 


8 


3 


.61±1 


.01 


3 


.59 


14 


.19 


+0 


' 54 -0 


30 

.55 


— 0. 


07 -0 


25 


cxou 


T 1 OO E O A 

,J IZZbZu 


. I + ls(J9bb 


12 


25 


20 


. 72 


18 


09 


55 


.7 


101. 


8 


4 


.65±1 


.14 


4 


.09 


18 


.29 


— 


12 -0 


39 
.36 


— 0. 


,.+0. 
35 -0 


26 


X u 


J122520 


. o+ loioU 1 


12 


25 


20 


.32 


18 


13 


01 
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Notes 

(ii) 



Note. - 



The notation is the same as in Table 1. 



Table 3. X-ray Spectral Fits of NGC 4365 



Hard Component Soft Component (mckal) 

N H Model kT h or F F% kT a Abund. Ff) 

Row Origin Region (10 2() cm" 2 ) (keV) ( a ) (kcV) (solar) ( a ) Counts X 2 /d°f 
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Sources 


< 1 O-eff 


(1.63) 


Brcmss 


6.08_! 62 
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Sources 


< 1 a cB 


U.UU [<.O.DOJ 
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„ --+3.11 
b - 5 '-1.83 


qb 
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< 1 left 


(1.63) 


Power 
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1.67_ 12 
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u.uu [^y. 1 4j 
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iu.o*_ 10 2g 


Power 


-1 7 o+0.24 
1 - lo -0.22 
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Field 
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(1.67) 
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Field 


(1.63) 


Power 


1 -°°-0 .06 
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Sources 


Field 


3.04 [<8.76] 


Power 
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Unresolved 
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Unresolved 


< 1 a e s 


0.00 [<13.11] 






13 b,c 


Unresolved 


< 1 a cB 
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887 


oo c\ /on "i -to 
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887 
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a Units arc 10 13 ergs cm 2 s 1 in 0.3-10 keV band. 

b The adopted best-fit model for this emission. 

c The energy range for this spectra excludes 1.6 - 1.9 keV. 



Table 4. X-ray Spectral Fits of NGC 4382 
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Note. - 



Notation is the same as Table 3. 
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Fig. 1.— Adaptively smoothed Chandra S3 image (0.3-10 keV) of (left) NGC 4365 and 
(right) NGC 4382, corrected for exposure and background. The gray scale is logarithmic 
and ranges from 1.56 x 10~ 6 cnt pix -1 s _1 to 3.3 x 10~ 4 cnt pix -1 s _1 . 




Fig. 2.— Logarithmic gray scale DSS optical image of (left) NGC 4365 and (right) NGC 4382. 
The circles indicate the positions of the X-ray sources listed in Tables 1 and 2. The relative 
surface brightness range covered by the gray scale is narrower than in Figure 1. 
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Fig. 3. — Histogram of the observed cumulative luminosity function of all resolved sources 
within the Chandra S3 fields of (left) NGC 4365 and (right) NGC 4382. The continuous 
curves are the sums of the respective best-fit LMXB luminosity functions (eq. 1) and the 
expected background source counts. Poisson error bars are displayed for every fifth interval. 
As this is a cumulative distribution, the errors are correlated. 
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Fig. 4. — Hardness ratios for the sources of (left) NGC 4365 and (right) NGC 4382 with at 
least 20 net counts. Here, H21 = (M — S)/(M + S) and #31 = (H — S)/(H + S), where S, 
M, and H are the counts in the soft (0.3-1 keV), medium (1-2 keV), and hard (2-10 keV) 
bands, respectively. The area of each circle is proportional to the observed number of net 
counts. The solid curve and large diamonds show the hardness ratios for power-law spectral 
models with Galactic absorption; the small diamonds show the ratios for N H = 5 x 10 21 
cm~ 2 ; the diamonds indicate values of the power-law photon number index of T = (upper 
right) to 3.2 (lower left) in increments of 0.4. The dashed curve encloses the area with 
power-law photon indices between and 3.2 and absorption columns between Galactic and 
5 x 10 21 cm" 2 . The la error bars at the upper left illustrate the median of the uncertainties. 
Due to the QE degradation, an effective absorption of ~ 6 x 10 20 cm" 2 was added to these 
spectral models. 
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Fig. 5. — Solid histogram gives accumulated fraction of events for sources as a function 
of exposure time. The dashed line is the predicted distribution under the hypothesis that 
the source plus background rate is constant. The three NGC 4365 sources shown are those 
with variability detected at the > 99% confidence level. NGC 4382 Src. 54 had variability 
detected at the > 96.8% confidence level. 
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Fig. 6. — The X-ray light curve of Src. 86 in NGC 4365. The points with la error bars are 
the source count rate, uncorrected for background, accumulated in 2021 s time bins. The 
horizontal dashed line is the background count rate scaled to the source region size. 
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Fig. 7. — Solid histogram gives the distribution of the position angles of the X-ray sources 
within 3' of the center of (left) NGC 4365 and (right) NGC 4382 in 10° bins, as a function 
of the PA (modulo 180°). PA is measured from north to east. The curves are the expected 
number of background sources plus the predicted distribution based on optical photometry 
and the best-fit elliptical distribution determined by a maximum likelihood fit to the observed 
values. 
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Fig. 8. — Solid histogram gives the accumulated number of X-ray sources within a 3' el- 
liptical isophote of the center of (left) NGC 4365 and (right) NGC 4382 as a function of 
semi-major axis. The curves are the expected number of background sources plus the pre- 
dicted distribution based on optical photometry and the best-fit distribution determined by 
a maximum likelihood fit to the observed values. Poisson error bars are displayed for every 
fifth interval. As this is a cumulative distribution, the la error bars are correlated. 
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Fig. 9. — In the upper panels, the cumulative X-ray spectra of all of the resolved sources 
in the inner effective radius of NGC 4365 (left) and inner two effective radii of NGC 4382 
(right) with la error bars are overlayed by the solid histograms of the best-fit model spectra 
(Table 3 and 4 row 3). Lower panels display the contribution to x 2 with the sign indicating 
the sign of the residual. 
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Fig. 10. — In the upper panels, the cumulative X-ray spectra of the unresolved emission 
in the inner effective radius of NGC 4365 (left) and inner two effective radii of NGC 4382 
(right) with la error bars are overlayed by the solid histograms of the best-fit model spectra 
(Table 3 row 13 and 4 row 11). Lower panels display the contribution to x 2 with the sign 
indicating the sign of the residual. 
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Fig. 11. — Surface brightness profiles of the soft band (0.3-1 keV) diffuse gaseous emission of 
(left) NGC 4365 and (right) NGC 4382 as a function of projected radius r with la error bars. 
The dashed curve shows the RC3 de Vaucouleurs profile with an effective radius determined 
from the distribution of optical light in the galaxy, but with the normalization varied to fit 
the X-ray surface brightness. The solid curve is the best-fit beta model. All fits were for 
r < 2'. 



